The muon g-2 and the bounds on the Higgs boson mass 
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After a brief review of the muon g— 2 status, we analyze the possibility that the present discrepancy 
between experiment and the Standard Model (SM) prediction may be due to hypothetical errors in 
the determination of the hadronic leading-order contribution to the latter. In particular, we show 
how an increase of the hadro-production cross section in low-energy e + e~ collisions could bridge the 
muon g— 2 discrepancy, leading however to a decrease on the electroweak upper bound on M H , the 
SM Higgs boson mass. That bound is currently M H 150 GeV (95%CL) based on the preliminary 

top quark mass M t = 172.6(1.4) GeV and the recent determination Ao^(M z ) = 0.02768(22), while 
the direct-search lower bound is M H > 114.4 GeV (95%CL). By means of a detailed analysis we 
conclude that this solution of the muon g— 2 discrepancy is unlikely in view of current experimental 
error estimates. However, if this turns out to be the solution, the 95%CL upper bound on M H 
is reduced to about 130 GeV which, in conjunction with the experimental lower bound, leaves a 
narrow window for the mass of this fundamental particle. 

PACS numbers: 13.40.Em, 14.60.Ef, 12.15.Lk, 14.80.Bn 



A. Introduction 

The measurement of the anomalous magnetic moment 
of the muon a M by the E821 experiment at Brookhaven, 
with a remarkable relative precision of 0.5 parts per mil- 
lion [1], is challenging the Standard Model (SM) of par- 
ticle physics. Indeed, as each sector of the SM con- 
tributes in a significant way to the theoretical prediction 
of a M = (g — 2)/2 (g is the muon's gyromagnetic factor), 
this measurement allows us to test the entire SM and pro- 
vides a powerful tool to scrutinize viable "new physics" 
appendages to this theory [2], 

The SM prediction of the muon g— 2 is conveniently 
split into QED, electroweak (EW) and hadronic (leading- 
and higher-order) contributions: a„ M = a^ ED + a EW + 
fl HLo _|_ a HHo rp^ e ^(Jj-qj^q contributions dominate the 
present a™ uncertainty. The QED prediction, com- 
puted up to four (and estimated at five) loops, currently 
stands at a« ED = 116584718.10(16) x 10~ n [3, 4], while 
the EW effects, suppressed by a factor (m M /M w ) 2 , pro- 
vide a EW = 154(2) x lO^ 11 [5]. The most recent calcula- 
tions of the hadronic leading-order contribution via the 
hadronic e + e~ annihilation data, to be discussed later, 
are in very good agreement: a" LO = 6909(44) x 10 _11 [6], 
6894(46) x 10- n [7], 6921(56) x l(T n [8], and 6944(49) x 
10 _11 [9]. The higher-order hadronic term is further di- 
vided into two parts: a™° = a™°(vp) + a™°(lbl). The 
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first one, -98(1) x 10 _11 [7], is the 0{a 3 ) contribution 
of diagrams containing hadronic vacuum polarization in- 
sertions [10]. The second term, also of 0(a 3 ), is the 
hadronic light-by-light contribution; as it cannot be de- 
termined from data, its evaluation relies on specific mod- 
els. Recent determinations of this term vary between 
80(40) x 10- n [ll] and 136(25) x 10~ n [12]. The most 
recent one, 110(40) x 10 _11 [13], lies between them. If 
we add this result to the leading-order hadronic contri- 
bution, for example the value of Ref. [7] (which also pro- 
vides a recent calculation of the hadronic contribution 
to the effective fine-structure constant, later required for 
our analysis), and the rest of the SM contributions, we 
obtain a" M = 116591778(61) x 10" 11 . The difference with 
the experimental value a EXP = 116592080(63) x lO" 11 [1] 
is Aa M = ajf p - al M = +302(88) x HT 11 , i.e., 3.4 stan- 
dard deviations (all errors were added in quadrature). 
Similar discrepancies are obtained employing the values 
of the leading-order hadronic contribution reported in 
Refs. [6, 8, 9]. 

The term a™° can alternatively be computed in- 
corporating hadronic r-decay data, related to those of 
hadroproduction in e + e~ collisions via isospin symme- 
try [14, 15]. Unfortunately there is a large difference 
between the e + e~- and r-based determinations of a" LO , 
even if isospin violation corrections are taken into ac- 
count [16]. The r-based value is significantly higher, 
leading to a small (~ la) Aa M difference. As the e + e~ 
data are more directly related to the a™° calculation 
than the r ones, the latest analyses do not include the 
latter. Also, we note that recently studied additional 
isospin-breaking corrections somewhat reduce the differ- 
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ence between these two sets of data (lowering the r-based 
determination) [17, 18], and a new analysis of the pion 
form factor claims that the r and e + e~ data are con- 
sistent after isospin violation effects and vector meson 
mixings are considered [19]. Recent reviews of the muon 
5 -2 can be found in Refs. [4, 20, 21]. 

The 3.4 a discrepancy between the theoretical predic- 
tion and the experimental value of the muon g—2 can 
be explained in several ways. It could be due, at least 
in part, to an error in the determination of the hadronic 
light-by-light contribution. However, if this were the only 
cause of the discrepancy, a" HO (lbl) would have to move 
up by many standard deviations to fix it - roughly eight, 
if we use the a™°(lbl) result of Ref. [13] (which includes 
all known uncertainties), and more than ten if the less 
conservative estimate of Rcf. [12] is employed instead. 
Although the errors assigned to a™°(lbl) are only edu- 
cated guesses, this solution seems unlikely, at least as the 
dominant one. 

Another possibility is to explain the discrepancy Aa M 
via the QED, EW and hadronic higher-order vacuum po- 
larization contributions; this looks very improbable, as 
one can immediately conclude inspecting their values 
and uncertainties reported above. If we assume that the 
g— 2 experiment E821 is correct, we are left with two op- 
tions: possible contributions of physics beyond the SM, or 
an erroneous determination of the leading-order hadronic 
contribution a^ LO (or combinations of the two). The first 
of these two options has been widely discussed in the lit- 
erature; we will focus on the second one, and analyze its 
implications for the EW bounds on the mass of the Higgs 
boson. 



B. Shifts of a£ LO and Aq^(M z ) 

The evaluation of the hadronic leading-order contri- 
bution a^ LO , due to the hadronic vacuum polarization 
correction to the one-loop QED diagram, involves long- 
distance QCD for which perturbation theory cannot be 
employed. However, using analyticity and unitarity, it 
was shown long ago that this term can be computed from 
hadronic e + e~ annihilation data via the dispersion inte- 
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where a(s) is the total cross section for e + e~ annihilation 
into any hadronic state, with extraneous QED corrections 
subtracted off, and s is the squared momentum transfer. 
The kernel K(s) is the well-known function 
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monotonically for increasing s and, for large s, it be- 
haves as rn 2 ,/ (3s) to a good approximation. One finds 
that the low-energy region of the dispersive integral is 
enhanced by ~ 1/s 2 . About 90% of the total contribu- 
tion to a™° is accumulated at center-of-mass energies 
yfs below 1.8 GeV and roughly three- fourths of a^ LO is 
covered by the two-pion final state which is dominated 
by the p(770) resonance [15]. Note that a™° is a pos- 
itive definite quantity. Exclusive low-energy e + e~ cross 
sections have been measured by experiments running at 
e+e~ colliders in Frascati, Novosibirsk, Orsay, and Stan- 
ford, while at higher energies the total cross section has 
been measured inclusively. Perturbative QCD becomes 
applicable at higher loop-momenta, so that at some en- 
ergy scale one can switch from data to QCD [24]. 

Let's now assume that the discrepancy Aa M = a® XF — 
a™ = +302(88) x 10" 11 , is due to - and only to - hy- 
pothetical mistakes in cr(s), and let us increase this cross 
section in order to raise a" LO , thus reducing Aa p . This 
simple assumption leads to interesting consequences. An 
upward shift of the hadronic cross section also induces 
an increase of the value of the hadronic contribution to 
the effective fine-structure constant at Mz [25], 
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(P stands for Cauchy's principal value). This integral is 
similar to the one we encountered in Eq. (1) for a" L °. 
There, however, the weight function in the integrand 
gives a stronger weight to low-energy data. The negligible 
contribution to a" LO and Aa^ d (M z ) of the 7r°7 channel 
below the n + n~ threshold was ignored in Eqs. (1,3). Let 
us define 



a = dsf(s)a(s), 

Ural 



dsg{s)cr(s), 



(4) 
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where the upper limit of integration is s u < M 2 , and the 
kernels are f(s) = K(s)/(4tt 3 ) and g(s) = [M 2 /(M 2 - 
s)]/(4a7r 2 ). Equations (4,5) provide the contributions to 
a™° and Aa^ d (M z ), respectively, in the region from the 
two-pion threshold up to s u (see Eqs. (1,3)). 

An increase of the cross section a(s) of the form 



Aer(s) = ecr(s) 



(6) 



in the energy range y/s e [\/s — 5/2, y/s +S/2], where e is 
a positive constant and 2m^+S/2 < t/s < \fs u — S/2, in- 



creases a by Aa(y / s , S, e) = e 



\/5 -«/2 



2ta(t 2 )f(t 2 )dt. 

If we assume that the muon g—2 discrepancy is entirely 
due to this increase in <r(s) so that Aa(^/s , 5, e) = Aa^, 
the parameter e becomes 



(see Ref. [23] for some of its explicit representations and 
their suitability for numerical evaluations). It decreases 



(7) 
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and the corresponding increase in Aa^ d (M z ) is 

A6(Vi , S) = Aa /^~ + 5 , (8) 

In the limiting case of a point-like shift Aa(s) = 
e'S(s — so), with 2m ff < v^o ^ V^ui the condition 
Aa(y/s ,e') = Aa M , with Aa( v / i ,e') = e'f(s ), leads 
to 

A6(^ ) = Aa„ [9(« )//(*o)] • (9) 

Following Ref. [15], to overcome the lack of precise data 
for a(s) at threshold energies, in the region 2m w < ^/s < 
500 MeV one can adopt the polynomial parametrization 
for the pion form factor F n (s) inspired by chiral pertur- 
bation theory; the parameters are determined from a fit 
to the data for both time-like and space-like momentum 
transfers [15, 23, 26]. The cross section below 500 MeV 
is therefore given by 

a(s) = ?—ft\F 7r (s)\ 2 , (10) 

6S 

where /?„ = (1 - 4m 2 /s) 1 / 2 , F^s) = 1 + s(r 2 )^/6 + 
s 2 ci + s 3 c 2 , (r 2 )^ = (0.439 ± 0.008) fm 2 , ci = (6.8 ± 
1.9) GeV~ 4 , and c 2 = (-0.7 ± 6.8) GeV~ 6 (see Ref. [15] 
for the correlation matrix of these coefficients). Be- 
tween 500 MeV and 1.4 GeV we use the cross section di- 
rectly obtained combining the experimental results of the 
tt+tt 1 [27], tt+tt-tt [28, 29], K + K~ [29, 30], K° L K° S [31], 
27T+27T- [32], 7r°7r°7r+7r- [33], 7r° 7 [34, 35] and 777 [35] 
channels. Between 1.4 GeV and 2 GeV we employ the 
inclusive measurements of Ref s. [36]. 

Figure 1 shows the shifts Ab(-y/s ,S = 210MeV) (his- 
togram) and A6(y / s ) (smooth curve) obtained from the 
increases Acx(s) = ecr(s) and Aa(s) — e'S(s — s ), respec- 
tively. These shifts, shown as functions of v^o' are a dded 
to the value Aa^ d (M z ) = 0.02768(22) [7]. The uncer- 
tainty of the sum Aa^ d (M z ) + Ab(^/s Q ) is indicated by 
the light band. To compute it, we first note that the er- 
rors 46x lO" 11 in and 22x 10~ 5 in Aa^(M z ) [7] are 
strongly correlated since they arise mainly from the same 
source, namely the uncertainty in the hadronic e + e~ an- 
nihilation cross section (which includes the uncertainties 
associated with the radiative corrections applied to the 
experimental data). Taking this into account, and ob- 
serving also that the error in Ab(y/s Q ) due to the aJJ LO 
uncertainty is —46 x 10~ n [<?(so)//(so)]j we add it lin- 
early to 22 x 10~ 5 , and then combine in quadrature this 
result with the error in A6(y / s ) induced by the remain- 
ing Aa M uncertainty. (We note that combining all errors 
in quadrature, ignoring their correlation, would enlarge 
the uncertainty of the sum Aaj^ d (M z ) + A6(y / s ), but 
would only induce minimal changes in our analysis.) The 

uncertainty of the sum Aa^ d (Ad z ) + Ab(^/s Q: 5), for fi- 
nite energy intervals, is computed analogously, neglecting 
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FIG. 1: Shifts of Aa^(M z ). The histogram indicates the 
increase Ab(i/s ,S) obtained varying the cross section by 
Acr(s) = eer(s) in 5 = 210 MeV energy regions, while A6(y / s ), 
obtained for point-like increases, is plotted as a smooth curve. 
The shifts are added to Aa^ d (M z ) = 0.02768(22) [7] (hor- 
izontal line). The uncertainty of the sum Aa[^ d (M z ) + 
A6(y / s ) is shown by the light band. 



the relative error of the ratio of integrals on the r.h.s. of 
Eq. (8) with respect to the large relative error of Aa M . 
The dark area below 2to„., where is the mass of the 
charged pion, denotes the kinematically forbidden region 
below the ir + ir~ threshold (the ir°j channel is neglected 
below this threshold). 



C. Connection with the Higgs boson mass 

The dependence of SM predictions, via quantum ef- 
fects, on the mass of the Higgs boson M H provides a 
powerful tool to set indirect bounds on the mass of 
this fundamental missing piece of the SM. Indeed, com- 
paring calculated quantities with their precise experi- 
mental values, the present global fit of the LEP Elec- 
troweak Working Group (LEP-EWWG) leads to the value 
M H = 87la7 GeV and to the 95% confidence level 
(CL) upper bound Mf ~ 160 GeV [37]. This result 
is based on the very recent preliminary top quark mass 
M t = 172.6(1.4) GeV from a combined CDF-DO fit [38] 

and the value Aa[^ d (M z ) = 0.02758(35) [39]. The LEP 
direct-search lower bound is M^ B = 114.4 GeV [40], also 
at the 95% CL. 

Although the global fit to the EW data employs a large 
set of observables, the M H upper bound is strongly driven 
by the comparison of the theoretical predictions of the 
mass of the W boson and the effective EW mixing an- 
gle sin 2 #^ t with their precisely measured values [41]. 
Convenient formulae providing the SM prediction of M w 
and sin^'g 5 * in terms of M H , the top quark mass M t , 
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Aa^ d (M z ), and a s (M z ), the value of the strong cou- 
pling constant at the scale M z , are given in Refs. [42]. 
Combining these M w and sin 2 0^ * predictions by means 
of a numerical x 2 -analysis, and using the present world- 
average values M w = 80.398(25) GeV [43-45], sin 2 ^ e ff pt = 
0.23153(16) [46], M t = 172.6(1.4) GeV [38], a s {M z ) = 
0.118(2) [47], and the determination Aa[ 5 a } d (M z ) = 
0.02758(35) [39] adopted by the LEP-EWWG, we obtain 
M H = 92±|| GeV and M 9 J' = 161 GeV. We see that in- 
deed the M H values obtained from the M w and sin 2 ^ 1 ^' 
predictions are quite close to the results of the global 
analysis. 

The M H dependence of the SM prediction of the muon 
g— 2, via its EW contribution, is too weak to provide 
M H bounds from the comparison with the measured 
value. Indeed, the shift of a™ for M H varying be- 
tween 114.4 GeV and 300 GeV is only of O(10- n ), 
which is negligible when compared with the hadronic 
and experimental uncertainties. On the other hand, 
Aa[^ d (M z ) is one of the key inputs of the EW fits. 
For example, employing the recent (slightly higher) value 
Aa^ d (M z ) = 0.02768(22) [7] instead of Aa ( ^ d (M z ) = 
0.02758(35) [39], the M H prediction shifts down to M H = 
90llg GeV and M h = 150 GeV. We note that M™ de- 
pends both on the central value and on the uncertainty 

of Aaj^ d (M z ). Henceforth, we employ the recent evalua- 
tion Aa [ ^ d (M z ) = 0.02768(22) [7]. (For the dependence 
of M H and its bounds on Aa^ d (M z ) see Refs. [42]). Next 

we consider the new values of Aa^ d (M z ) obtained shift- 
ing 0.02768(22) by Ab(^ ) and A6(V? , S) (including 
their uncertainties as discussed in the previous section), 
and compute the corresponding new values of M 9 ^ by 
means of the combined x 2_anar y s i s based on the M w 
and sin^gg?' inputs. The results are shown in Fig. 2. 
The lower region M H < 114.4 GeV is excluded by the di- 
rect LEP searches at 95% CL, while the upper one is ex- 
cluded by the indirect EW 95% CL bound M H < 150 GeV 
obtained with Aa { ^ d (M z ) = 0.02768(22). (As in the 

case of Aa^ d (M z ), the value adopted here for a™° 
is from the recent article in Ref. [7].) If we increase 
the hadronic cross section a(s) by e'S(s — so) in or- 
der to bridge the muon g— 2 discrepancy Aa^, M^ 5 de- 
creases, as shown by the continuous red line in Fig. 
2, further restricting the already narrow allowed region 
for M H . In particular, this curve falls below M^ B for 
\/s ^1.1 GeV. The two histograms show the M^ 5 val- 
ues when the analysis is repeated with Act = ecr(s) shifts 
in 5 = 210 MeV and 5 = 400 MeV energy regions. We 
conclude that the hypothetical shifts Act = eer(s) (in 
s/s 6 [-y/s — S/2, y / s + (5/2]) of the hadronic cross section 
that bridge the muon g— 2 discrepancy, conflict with the 
LEP lower limit when y/s Q > (\/s )thr ~ 1.2 GeV, for 
values of S up to several hundreds of MeV. The thresh- 
old (\/so)thr increases above ~ 1.3 GeV for hypothetical 
shifts ecr(s) in even wider energy regions S 1 GeV, but 
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FIG. 2: The Mf values obtained via the M w and sin 2 0*ff 
fits using as input for Aa^(M z ) the value 0.02768(22) in- 
creased by A6( v / s ) (smooth curve) and by Ab(y r s Q ,5 = 
210MeV,400MeV) (histograms). The area below 114.4 GeV, 
partly yellow and partly pink, is excluded at 95% CL by the 
LEP direct lower bound, while the orange M H > 150 GeV 
one is forbidden by the EW indirect upper bound. As in 
Fig. 1, the region y/s < 2m^ is excluded. The dotted line 
replaces the smooth one when r data are incorporated in the 
determination of Aa^ d (M z ) and a^ M . 

uniform shifts of the cross section in such wide energy 
ranges appear to be unrealistic. 

If t data are incorporated in the calculation of the 
dispersive integrals in Eqs. (1) and (3), the leading- 
order hadronic contribution to the muon g— 2 signifi- 
cantly increases to a« LO = 7110(58) x 10~ n [15], the 
higher-order vacuum polarization term slightly decreases 
to a™°(vp) = -101(1) x 10~ n [7, 20], and the dis- 
crepancy with the experimental value drops to Aa M = 
+89(95) x 10~ n , i.e. roughly 1 a. While using r data 
almost solves the muon g— 2 discrepancy, it increases the 
value of Aa^{M z ) to 0.02782(16) [15, 48]. In Ref. [48] it 
was shown that this increase leads to a low M H prediction 
which is suggestive of a near conflict with M^ B , leaving 
a very narrow window for M H . Indeed, with this value of 

Aa hld( M z) and thc 

same above-discussed values of the 
other inputs of the x 2 -analysis, we find M H = 84^3 GeV 
and an M^ 5 value of only 138 GeV. The dotted line in 
Fig. 2 shows the M 9 ^ values obtained using r data to 

compute Aa^\{M z ) and Aa M , with the hadronic cross 
section a(s) increased by e'S(s — sq) in order to bridge 
the Aa M difference. 

As we briefly mentioned in the Introduction, recently 
computed isospin-breaking violations, further improve- 
ments of the long-distance radiative corrections to the 
decay t~ — ► ii~ti {) v t [17] and differentiation of the neu- 
tral and charged p properties [18], reduce to some extent 
the difference between r and e + e~ data, lowering the r- 
based determination of a" L °. Moreover, a recent analysis 
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of the pion form factor below 1 GeV claims that r data 
are consistent with the e + e~ ones after isospin violation 
effects and vector meson mixings are considered [19]. In 
this case one could therefore use the e + e~ data below 
~ 1 GeV, confirmed by the r ones, and assume that Aa^ 
is accommodated by hypothetical errors occurring above 
~ 1 GeV, where disagreement persists between these two 
sets of data. Our previous analysis shows that this as- 
sumption would lead to values of M^ 5 inconsistent with 
the LEP lower bound. 

It is interesting to note that there are more complex 
scenarios where it is possible to bridge the Aa^ discrep- 
ancy without significantly affecting M 9 J" . For instance, 
we may envisage an increase of <r(s) at low s combined 
with a decrease at high s in such a manner that their 

overall contribution to Aa^ d (M z ), and therefore to M^ 5 , 
approximately cancels. Since the contributions to a^ LO 
are more heavily weighted at low s, it is then possible 
to further adjust the positive and negative a(s) shifts to 
bridge the muon g— 2 discrepancy. However, such scheme 
requires two fine-tuning steps and a larger increase of 
ct(s) at low s, and is therefore considerably more unlikely 
than the simplest scenarios, involving a single adjustable 
contribution, that are discussed in detail in this paper. 



the shift Act = ecr(s) in that energy range. A shift up 
of a(s) in the energy range from 2771^ to 850 MeV, to 
fix Aa M , leads to e - 6% and lowers to 134 GeV. 
Higher values of e are obtained for narrower energy bins, 
particularly if they do not include the p-uj resonance re- 
gion. For example, a huge e ~ 52% increase is needed 
to accommodate Aa M with a shift of the cross section in 
the region from 2m 7r up to 500 MeV (reducing M^ 5 to 
143 GeV), while an increase in a bin of the same size but 
centered at the p peak requires e ~ 8% (lowering M^ 5 
to 132 GeV). As the quoted experimental uncertainty 
of <r(s) below 1 GeV is of the order of a few per cent 
(or less, in some specific energy regions), the possibility 
to explain the muon g— 2 discrepancy with these shifts 
Acr(s) appears to be unlikely. Figure 3 shows that for 
fixed 5 (i.e., segments of the same color), lower values of 
e are obtained if the shifts occur in energy ranges cen- 
tered around the p-to resonances; but also this possibility 
looks unlikely, since it requires variations of ct(s) of at 
least ~ 6%. If, however, we allow variations of the cross 
section up to ~ 6% (7%), is reduced to less than 
- 134 GeV (135 GeV). For example, the - 6% shifts in 
the intervals [0.5,1.0] GeV or [0.6,1.2] GeV, required to 
fix Ao^ (not represented in Fig. 3), lower to 133 GeV 
or 130 GeV, respectively. 



D. How realistic are these shifts A<r(s)? 

In the above study, the hadronic cross section a(s) 
was shifted up by amounts required to adjust the muon 
g— 2 discrepancy Aa M . Apart from the implications for 
the Higgs boson mass (and the restrictions deriving from 
them), these shifts may actually be inadmissibly large 
when compared with the quoted experimental uncertain- 
ties. For example, one of the histograms in Fig. 2 shows 
that a shift Act in a 210 MeV bin centered just above 
the p peak could fix the muon g— 2 discrepancy (lower- 
ing M^ 5 to 131 GeV); but is such a shift of the precisely 
measured cross section at the p peak realistic? 

To investigate this problem, we turn our attention to 
the parameter e = Act(s)/ct(s), i.e. the ratio of the shift 
Act(s) required to bridge the muon g— 2 discrepancy and 
the cross section ct(s), provided by Eq. (7). Clearly, the 
value of e depends on the choice of the energy range 
[y/s Q — (5/2, y/s + 5/2] where ct(s) is increased and, for 
fixed y/s , it increases when 5 decreases. The minimum 
value of e is roughly +4%; it occurs if the hadronic cross 
section ct(s) is multiplied by (1 + e) in the whole integra- 
tion region of Eq. (1), from the tt + it~ threshold to infinity 
(this minimum value of e changes only negligibly whether 
the shift up of ct(s) includes or not the high-energy re- 
gion where perturbative QCD is employed). Such a shift 
would lead to an M 9 ^ of roughly 75 GeV, well below the 
LEP lower bound. 

Figure 3 shows the values of e (in per cent) for several 
bin widths 5 and central values y/s (same length seg- 
ments are of the same color). Also, next to each segment 
we quote the value of M^ 5 (in GeV) obtained performing 
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FIG. 3: Values of e obtained increasing er(s) by ecr(s), 
to bridge the Aa^ discrepancy, in energy ranges [y^o ~~ 
5/2, \fs + 5/2] for various values of y/s and <5. The number 
next to each segment indicates the M'g value (in GeV) in- 
duced by the eer(s) shift in that energy region. Same length 
segments are of the same color. The midpoint of each segment 
is displayed by a dot. 

We remind the reader that the present experimental re- 
sults for sin 2 ^^' exhibit an intriguing dichotomy. Those 

based on the leptonic observables lead to (sin 2 ^ 1 ^*); = 
0.23113(21), while the average of those derived from the 
hadronic sector is (sin 2 ^° ff pt ) h = 0.23222(27) [46]. The 
results within each group agree well with each other, but 
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the averages of the two sectors differ by about 3.2a. Our 
analysis, like the LEP-EWWG one, depends on the value 
of sin 2 ^^?*. For instance, if we were to use (sin 2 ^^*)/,, 
we would obtain a significantly higher SM prediction: 
M H = 129±g GeV, Mf = 225 GeV, and a continu- 
ous (red) line in Fig. 2 similarly shifted up. However, we 
note that in this scenario the M H predictions from M w 
and (sin^'g 5 *)/! are inconsistent with one another unless 
one introduces additional "new physics" beyond the SM. 
For example, the difference could be associated with a 
value S <~ 0.4 to 0.5 of the S'-parameter, an effect gener- 
ally attributed to technicolor-like theories with additional 
heavy fermion chiral doublets [49] . Instead, if we were to 
employ (sin 2 ^^?*)/, the SM prediction would drop sharply 
to M H = 501 2 ] GeV, M 9 H 5 = 97 GeV, which is already 
in conflict with the direct-search lower bound. Thus, in 
that case, no shift Aa(s) could reconcile the g— 2 discrep- 
ancy without violating the lower bound. In this paper we 
employ as input the world-average of sin 2 6^* since this 
is the value determined in the global analysis of the SM. 

The M H upper bounds presented in this article depend 
sensitively on the central value M t = 172. 6 GeV and its 
uncertainty SM t = 1.4 GeV. In the future, the former 
may still change and the latter will further decrease. We 
therefore provide the following simple formulae to trans- 
late easily the M^ 5 = 150 GeV result of our numerical 
% 2 -analysis based on the M w and sin 2 ^^' predictions, 
as well as the M^ 5 [0.6, 1.2] = 130 GeV upper bound cor- 
responding to the <~ 6% increase of a(s) in the interval 
[0.6,1.2] GeV (an illustrative case that accounts for Aa M ), 
into the new values derived with different M t and SM t 
inputs: 

M 9 H 5 = (150.5+ 11. 2.x + 9.4y) GeV, (11) 
M° 5 [0.6,1.2] = (130.7 + 9.9a; + 8. 2y) GeV, (12) 

with x = M t - 172.6GcV and y = SM t - 1.4GeV. Note 
that, in case of a future rise of the M t central value, 
the increase induced on the M H upper bounds would be 
partially compensated by a reduction of the error 6M t . 
Equations (11,12) reproduce the results of the detailed 
numerical x 2 -analysis with maximum absolute deviations 
of roughly 1 GeV when M t e [171, 174] GeV and SM t € 
[1.0,1.8] GeV. 



E. Conclusions 

The present discrepancy between the SM prediction of 
the anomalous magnetic moment of the muon and its ex- 
perimental determination could be due to the contribu- 
tion of new, yet undiscovered, physics beyond the SM, or 
to errors in the determination of the hadronic contribu- 
tions. In this letter we considered the second hypothesis 
and, in particular, the possibility to accommodate the 
discrepancy Aa M = +302(88) x 10~ n (3.4 a) by changes 
in the hadronic cross section a(s) used to determine the 



leading hadronic contribution a™°. This option has im- 
portant consequences on M^ 5 , the 95% CL EW upper 
bound on the mass of the SM Higgs boson. 

We first analyzed the effects induced by these hypo- 
thetical changes Acr(s) on the value of Aa^ d (M z ), one 
of the key inputs of the EW fits with a strong influence 
on the SM M H predictions. The comparison of the the- 
oretical predictions of M w and the effective EW mixing 
angle sin 2 ^^* with their precisely measured values al- 
lowed us to determine, via a combined x 2 analysis, the 
variations of M 9 ^ induced by the shifts Aer(s). We con- 
cluded that if the hadronic cross section is shifted up in 
energy regions centered above ~ 1.2 GeV to bridge the 
muon g— 2 discrepancy, the Higgs mass upper bound be- 
comes inconsistent with the LEP lower limit. 

If r-decay data are incorporated in the calculation of 
a* M , the discrepancy Aa M drops to +89(95) x 10 -11 . 
While this almost solves the muon <?— 2 discrepancy, 
it raises the value of Aa[^ d (M z ) leading to M 9 ^ — 
138 GeV, increasing the tension with the LEP lower 
bound. One could also consider a scenario, suggested 
by recent studies, where the r data confirm the e + e~ 
ones below ~ 1 GeV, while a discrepancy between them 
persists at higher energies. If, in this case, Aa M is recon- 
ciled by hypothetical errors above ~ 1 GeV, where the 
data sets disagree, one also finds values of M 9 ^ incon- 
sistent with the 114.4 GeV lower bound. For example, 
if a{s) is shifted in the interval [1.0,1.8] GeV, we obtain 
M 9 H 5 = 108 GeV. 

We then questioned the plausibility of the variations 
Aa(s) =ea(s) required to fix Aa M . Their amounts clearly 
depend on the energy regions chosen for the change, but 
we showed that they are generally very large when com- 
pared with the actual experimental uncertainties. Given 
the small experimental uncertainty of a(s) below 1 GeV, 
the possibility to bridge the muon g— 2 discrepancy with 
shifts of the hadronic cross section appears to be un- 
likely. Smaller values of e (for fixed bin-widths 6) are 
needed when the shifts occur in energy regions centered 
around the p-u resonances; but also this possibility looks 
unlikely since it requires variations of a(s) of at least 
~ 6%, a large modification given current experimental 
error estimates. However, if this turns out to be the so- 
lution of the Aa M discrepancy, we conclude that M 9 ^ is 
reduced to roughly 130 GeV which, in conjunction with 
the 114.4 GeV lower bound, leaves a narrow window for 
the mass of this fundamental particle. Simple formulae 
were also provided to translate M H upper bounds derived 
in this paper into new values corresponding to M t and 
5M t inputs different from those employed here. 

If the Aa M discrepancy is real, it points to "new 
physics", like low-energy supersymmetry. In fact, an 
intriguing explanation of Aa p is provided by some su- 
persymmetric models, where it is reconciled by the addi- 
tional contributions of supersymmetric partners [2] and 
one expects M H < 135 GeV for the mass of the lightest 
scalar [50]. If, instead, the deviation is caused by an in- 
correct leading-order hadronic contribution, it leads to a 
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larger AaJ lad (M z ) and, correspondingly, to low values of 
M^ 5 , thus leaving a very narrow range for the SM Higgs 
boson mass. 
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